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ABSTRACT 

We present results from a campaign of multiple epoch echelle spectroscopy of relatively faint 
(V = 9.5 — 13.5 mag) red giants observed as potential astrometric grid stars for the Space 
Interferometry Mission (SIM PlanetQuest). Data are analyzed for 775 stars selected from the 
Grid Giant Star Survey spanning a wide range of effective temperatures ( T e ff ), gravities and 
metallicities. The spectra are used to determine these stellar parameters and to monitor radial 
velocity (RV) variability at the 100 m s _1 level. The degree of RV variation measured for 489 
stars observed two or more times is explored as a function of the inferred stellar parameters. 
The percentage of radial velocity unstable stars is found to be very high - about 2/3 of our 
sample. It is found that the fraction of RV-stable red giants (at the 100 nr s -1 level) is higher 
among stars with T e ff~ 4500 K, corresponding to the calibration-independent range of infrared 
colors 0.59 < (J — K s )g < 0.73. A higher percentage of RV-stable stars is found if the additional 
constraints of surface gravity and metallicity ranges 2.3 < log g < 3.2 and -0.5 < [Fe/H] < 
-0.1, respectively, are applied. Selection of stars based on only photometric values of effective 
temperature (4300 K < T e ff < 4700 AT) is a simple and effective way to increase the fraction 
of RV-stable stars. The optimal selection of RV-stable stars, especially in the case when the 
Washington photometry is unavailable, can rely effectively on 2MASS colors constraint 0.59 < 
(J — K s ) o < 0.73. These results have important ramifications for the use of giant stars as 
astrometric references for the SIM PlanetQuest. 

Subject headings: stars: abundances, fundamental parameters, oscillations, late-type — techniques: ra¬ 
dial velocities 
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1. Introduction 

The Grid Giant Star Survey (GGSS; Patter¬ 
son et al. 2001) is a partially-filled, all-sky sur¬ 
vey to identify giant stars from 9 < V < 17.5 
mag using the Washington M, T-i + DDObl filter 
photometric prescription outlined in Majewski et 
al. (2000). A primary motivation for the survey 
and a driver of its design is the selection of stars 
suitable for the Astrometric Grid of the Space In- 
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terferometry Mission (SIM PlanetQuest). These 
Grid stars, which serve as astrometric references 
against which the motions of SIM targets are mea¬ 
sured, must themselves be as astrometrically sta¬ 
ble as possible. Thus, they must be free of stellar 
or significant planetary companions as well as at¬ 
mospheric activity (spotting/flaring) that will in¬ 
duce photocenter wobbles at the several microarc- 
second level on of order a decade-long timescale 
(approximately the duration of SIM). G and K 
giants were selected by the SIM project as the pri¬ 
mary stellar constituent of the SIM Astrometric 
Grid because these stars are the most luminous, 
common stellar type found in all directions of the 
sky. High intrinsic luminosity places giant stars 
at greater distances for a given apparent magni¬ 
tude, and this increased distance decreases the an¬ 
gular scale of any fixed linear astrometric wobble, 
thus making these particular stars less likely to be 
problematical as references. The GGSS project 
has been undertaken with the specific aim of iden¬ 
tifying subsolar metallicity giant stars, which are 
intrinsically more luminous than solar metallicity 
giants, and also less likely to have planets (see Fis¬ 
cher & Valenti 2005, and references therein). 

Despite consideration of such effects for care¬ 
ful selection, all SIM Astrometric Grid candidates 
require pre-mission monitoring to assess their like¬ 
lihood of astrometric stability during the mis¬ 
sion. One method is to search for variability in 
other properties, such as brightness and radial ve¬ 
locity (RV). Indeed, we are presently involved, 
along with other groups, in a large campaign of 
eclrcllc high-resolution RV monitoring of Astro- 
metric Grid candidates. Fortunately, not all vari¬ 
ability in RV (or luminosity) necessarily translates 
to detrimental astrometric variability. However, 
some intrinsic photometric and/or RV variability 
is due to processes that, although benign in terms 
of photocenter wobble (e.g., radial pulsations), can 
mask otherwise detectable signatures of problem¬ 
atical sources of variability, like astrometric wob¬ 
bles due to planetary, brown dwarf or stellar com¬ 
panions. Therefore, because campaigns to monitor 
radial velocity variability to the precision needed 
for vetting SIM targets (i.e. ~ 100 m s -1 ) are 
both expensive and time consuming, it is useful 
to understand how intrinsic RV variability of gi¬ 
ant stars depends on other stellar properties (like 
temperature, metallicity, and surface gravity) to 


optimize efficient selection of giant stars for the 
Astrometric Grid. 

Although it has long been known that stars at 
the top of the RGB are both brightness and RV- 
variable (see e.g. Pryor et al. 1988; Cote et al. 
1996), there has been little systematic study of 
this variability as a function of stellar atmospheric 
parameters, especially at fainter absolute magni¬ 
tudes. Monitoring of radial velocities for bright 
red giants has indicated a high probability of their 
RV variations at the level of 100 m s -1 (see Hatzes 
& Cochran 1994, and references therein). Joris- 
sen et al. (1997) found that red giants with spec¬ 
tral types late G to early K are stable and giants 
with later spectral types are all variable. At the 
same time, study of Hipparcos red clump giants 
(Adelman 2001) reveals photometric stability in 
this evolutionary phase. Radial velocity jitter at 
large amplitude is expected for metal-poor, lumi¬ 
nous red giants (My < -1.4) (see Carney et al. 
2003). 

In their assessment of the suitability of giant 
stars as SIM Astrometric Grid members, Frink et 
al. (2001) studied a proxy sample of 86 nearby gi¬ 
ant stars selected from the Hipparcos catalogue. 
Their study included three to eight epochs of high 
precision (5-8 m s -1 ) echelle RV measurements on 
timescales from days to a year. Frink et al. found 
that most (73 of 84) of the giant stars they in¬ 
vestigated have a stable RV at the level of < 100 
m s -1 , and the peak of the distribution of mea¬ 
sured dispersions in velocity occurs at 20 m s -1 . 
However, Frink et al. do report that RV varia¬ 
tion is more probable for redder stars ((B — V) > 
1.2 in their sample). While the closest match to 
the survey presented here, the Frink et al. survey 
contains stars that are much brighter than the ex¬ 
pected SIM Grid sample, and, moreover, they did 
not probe variability as a function of metallicity, 
which is a key aspect of the stars being found in 
the GGSS. 

Our goal here is to carry out an initial RV- 
stability assessment of stars more like those ex¬ 
pected to fill the Astrometric Grid by focusing on 
stars taken directly from the GGSS itself. The 
hope is that RV stability may correlate with some 
intrinsic stellar property, such as effective temper¬ 
ature, surface gravity, or metallicity. Although 
the precision of our RV measurements is an or¬ 
der of magnitude less than that of Frink et al. 
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(2001), it is appropriate for finding RV wobbles 
at the level needed to identify astrometrically- 
detrimental grid candidates for the fainter, > 1 
kpc distant GGSS candidates (and, indeed, is the 
precision at which monitoring campaigns of SIM 
Grid stars are being conducted). Moreover, our 
sample is almost an order of magnitude larger than 
that in the Frink et al. survey. The observations 
studied and discussed here were obtained under an 
initial JPL-sponsored program in which the inter¬ 
nal RV accuracy was set at about 100 ms -1 . 

2. Observations and Data Reduction 
2.1. Sample Selection and Biases 

The GGSS finds giant stars by photometry in 
the Washington M, T 2 + DDObl filters accord¬ 
ing to the methods described in Majewski et al. 
(2000). The (M - T 2 ,M - DD051) two-color 
diagram (2CD) effectively distinguishes between 
late type dwarfs and giants based on the surface- 
gravity-sensitive Mgb+MgH feature near 5150 A 
(where the 130A -wide DD051 filter is centered). 
Photometry was obtained in 1302 evenly placed 
fields across the celestial sphere, each of area 0.5- 
1.0 deg 2 . The 2CD (i.e., the Mgb+MgH feature) 
is secondarily sensitive to metallicity and the po¬ 
sitions of stars in the 2CD can therefore be used 
to derive a crude estimate of the metallicity of the 
likely giant stars. Metal-poor giants are the most 
widely separated stars in the 2CD from the locus 
of dwarf stars, and therefore these are the easiest 
to identify with our methods. Additional details 
regarding the GGSS can be found in Patterson et 
al. (2001). 

The selection of stars for the present study is 
biased according to the same specific set of cri¬ 
teria used in the original agreement of the GGSS 
collaboration with the SIM Project regarding the 
selection of Astrometric Grid candidates from the 
GGSS. 1 Giant star candidates are identified from 


1 In the past year, the SIM Project has modified somewhat 
the criteria for the selection of stars for the Astrometric 
Grid, with a focus on stars with magnitudes more typically 
V ~ 9 — 11 mag and chosen from both the GGSS and the 
Tycho-2 catalogues in a ratio of approximately 1:3. The Ty¬ 
cho stars are selected without foreknowledge of metallicity, 
and so are expected to be typically near solar metallicity 
and at smaller distances than stars from the GGSS. In this 
respect, all of our results are still relevant, however, since 
the stars included in the present study span the range of 


the 2CD and assigned photometric metallicities 
therefrom. From the apparent magnitude, color 
and metallicity, an estimate of the absolute mag¬ 
nitude and a photometric parallax distance is as¬ 
signed (Rhee et al. 2001). For each field, all giant 
candidates with V < 13.5 mag are rank-ordered 
according to distance, with the most distant star 
assigned highest priority. The top four ranked 
stars in each field, yielding more than 4000 can¬ 
didates of 9.5 < V < 13.5 mag, were passed 
onto both medium and high resolution spectro¬ 
scopic observing campaigns (the latter forming the 
database explored here). Typically, though not 
exclusively, the most distant giant candidates in 
each field were also among those with the lowest 
photometric metallicities brighter than V = 13.5 
mag. Thus, our spectroscopic sample is biased 
toward more metal-poor and distant stars than 
would be found from a random selection of gi¬ 
ant stars in the same magnitude range. From the 
above described sample of more than four thou¬ 
sand stars, we explore here a random subsample 
of 775 stars in both celestial hemispheres, which 
constitutes those stars for which multiple high- 
resolution spectra have been obtained during the 
first two years of the GGSS follow-up program. 

2.2. Spectroscopy 

Spectroscopic observations of 434 stars in the 
northern subsample of GGSS candidates were con¬ 
ducted in the period from January 2001 to De¬ 
cember 2002. We made use of the 2.1m telescope 
at the McDonald Observatory and the Sandiford 
Cassegrain echelle spectrograph, which provides 
R = 55000 resolution. Th-Ar comparison spec¬ 
tra were taken right before every program star ob¬ 
servation. The quality of the data was monitored 
by observing one or two RV standard stars from 
Nidever et al. (2002) per night. The adopted setup 
of the spectrograph enabled us to cover the 5000- 
5900 A spectral range and achieve a signal to noise 
(hereafter S/N) level of order 20-40 in reasonable 
exposure times (10-30 minutes). 

The observed spectra were reduced from two- 
dimensional to “echelle” format in a standard way 
with the IRAF 2 software package. Corrections for 

properties expected for both GGSS and Tycho-2 stars. 

2 IRAF is distributed by NOAO which is operated by AURA, 
Inc. under contract with the NSF. 
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bias, flat field, and scattered light were applied, 
and cosmic ray hits cleaned out by the IRAF’s 
tasks from crutil package. The internal accuracy 
of the wavelength calibration via the Th-Ar lamp 
spectra was on average 0.001 A (and not worse 
than 0.0018 A), which corresponds to an RV ac¬ 
curacy of order 55 m s^ 1 (and not worse than 100 
m s -1 ). 

Resampling the spectra introduces systematic 
errors in the wavelength calibration that degrades 
the RV accuracy. Thus, to preserve the wave¬ 
length calibration the spectrum orders were not 
concatenated to produce a single spectrum, but 
rather were maintained individually at their nat¬ 
ural sampling. 

Southern hemisphere GGSS stars were observed 
with the 1.2-rn Swiss telescope and CORALIE ve- 
locimeter in 2001-2004 at the ESO La Silla Ob¬ 
servatory (Chile), as described by Arenas et al. 
(2002). CORALIE is an improved version of the 
ELODIE spectrograph (Baranne et al. 1996). The 
effective resolution of CORALIE is 50000, it covers 
the wavelength range 3870 6800 A and it pro¬ 

vides a precision of typically 30 m s^ 1 in the radial 
velocities. Typically a S/N of 10 was achieved for 
program targets. The radial velocity and its ac¬ 
curacy comes directly from CORALIE’s reduction 
package TACOS (Baranne et al. 1996). 

CORALIE provides RVS for 375 program stars 
directly from the TACOS software. We found that 
spectra S/N < 4 have systematically low RV accu¬ 
racy and rejected them from further analysis. The 
whole sample of good S/N CORALIE candidates 
is thus reduced to 341 objects. 

While we cannot obtain abundances directly 
from CORALIE spectra, we incorporate the RVs 
of these stars in our analysis and adopt photomet¬ 
rically estimated stellar parameters for these stars. 

Since the southern sample was observed with a 
smaller aperture telescope than the northern one, 
it is biased toward brighter stars: the distribution 
of the southern V-magnitudes peaks at 11.6 mag 
whereas the northern sample distribution peak is 
located at 12.7 mag. 


3. Data Analysis 
3.1. Radial Velocities 

RVs from the McDonald spectra were measured 
using a cross-correlation methodology with the 
help of IRAF’s fxcor task. After the primary 
reduction the continuum of each order looks al¬ 
most flat except at its edges where the continuum 
level looks systematically lower and the S/N is de¬ 
graded. For the RV analysis we discard the pixels 
at order edges, which is about ~ 20% of the total. 

Each spectral order was cross-correlated against 
the RV template separately. The template cho¬ 
sen was the Arcturus spectrum (Hinkle et al. 
2000). The cross-correlated pieces of the tem¬ 
plate spectrum were convolved with a Gaussian 
corresponding to the FWHM of the Th-Ar lines 
found for each order. In addition, the corre¬ 
sponding piece of the calibrating Th-Ar spectrum 
was cross-correlated with a laboratory Th-Ar tem¬ 
plate (Palmer & Engleman 1983). The Th-Ar 
template for each order was also convolved with 
the same Gaussian as for the stellar template in 
that order. The shift in the radial velocity be¬ 
tween the observed Th-Ar and the laboratory ref¬ 
erence frame was subtracted from the stellar RV 
obtained from the cross-correlation via the tem¬ 
plate. Thus we not only have taken into account 
the varying instrumental profile by this order-by¬ 
order treatment, but we are able to make multiple 
estimates of the RV for a single observation and 
thereby assess the errors in our velocities for each 
star. A few of the spectral orders indicated sig¬ 
nificantly different RVs from the others because 
of poor wavelength calibration due to a sparse 
Th-Ar line sample. We discarded such ’’bad” or¬ 
ders from our analysis. The remaining orders’ ra¬ 
dial velocities ( RVi ) give us an independent as¬ 
sessment of the mean RV ( RVq ) and its accuracy: 
RV o = E RVi ■ Wi / E Wi, where Wi is the in¬ 
verse square of the individual RV error provided 
by IRAF. The dispersion of this estimate 5RV was 
found as 6RV = i/E w,. • (RVi — RVq) 2 / E Wi. 
We show an example of the order-to-order radial 
velocity difference assessed for the program star 
G1113+00.20 in Figure 1. 
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3.2. Stellar Parameters Derived from an 
Automated Spectroscopic Analysis 

In addition to the measured RV’s, physical stel¬ 
lar parameters can be derived for the red giants 
studied here. The effective temperature (T e ff) of 
giant stars can be determined from their broad¬ 
band colors. According to an observational cali¬ 
bration given by Alonso et al. (1999), near-infrared 
colors (J — K s ) can be used to estimate T e ff for 
red giants, even with no information on the stellar 
metallicity and surface gravity. The full release of 
the All-Sky Point Source Catalog of the Two Mi¬ 
cron All-Sky Survey (2MASS) makes available the 
near-infrared colors of all GGSS candidates. Fig¬ 
ure 2 shows the T e // distribution of the combined 
northern plus southern samples as estimated from 
(J — K s ) colors according to the Alonso et al. pre¬ 
scription. As may be seen, our candidates span a 
range of temperatures characteristic of red giants. 
The typical accuracy of 2MASS colors is 0.03 mag 
in our sample magnitude range. It corresponds to 
a 100 K accuracy in T e ff. 

The internal accuracy of the T e f f calibration in 
Alonso et al. (1999) from (J — K s ) is of order 125 
K. On the other hand, Bessell, Castelli, & Plez 
(1998) provide a calibration of the T e //(J — K s ) 
relation for giant stars from NMARCS stellar at¬ 
mosphere models. The Bessell, Castelli, & Plez 
results are reproduced by the Alonso et al. data 
fairly well for typical values of the surface grav¬ 
ity (see Figure 12 in Alonso et al. 1999). The 
relation from Bessell, Castelli, & Plez (1998) is 
metallicity-dependent but the difference between 
stars of [Fe/H]=0.0 and -2.0 (which corresponds 
to the bulk of our objects, see below) affects T e // 
by less than 100 K and would not have a signif¬ 
icant effect on the applied calibration procedure. 
From the previous two paragraphs, we expect ~ 
160 K accuracy in our derived T e f f. 

For the northern sample (434 stars) observed at 
McDonald Observatory, the reduced echelle spec¬ 
tra can be used to determine the surface grav¬ 
ity, log g, and metallicity, [Fe/H], by comparison 
to a library of synthetic stellar spectra. For this 
purpose we produced artificial spectra using the 
MOOG-2002 spectral synthesis package (Sneden 
1973). The corresponding stellar models were 
computed by the ATLAS 9 code (Kurucz 1993). 
We computed a model set covering the 3800 - 


5600 K temperature range, 0.2 < log g < 5.6, and 
—3.0 < [Fe/H] < +0.5 with steps in these param¬ 
eters of 200 K, 0.2 dex and 0.2 dex respectively. 
In all the models we assumed a typical value for 
the microturbulent velocity of 2 km s -1 . 

We selected three pieces of spectra to estimate 
the stellar parameters. The Arcturus spectrum 
(used as an RV template in §3.1) was split into 
40 A pieces, and the optimal synthetic spectrum’s 
piece and corresponding model parameters were 
found by least-square fitting. The stellar param¬ 
eters of Arcturus (obtained by Griffin & Lynas- 
Gray (1999)) were best reproduced in the regions 
of 5175 - 5215, 5215 - 5255, and 5700 - 5740 A. 
These wavelength intervals used to derive the stel¬ 
lar parameters are not simply dominated by Fe 
absorption lines, but also by (principally) lines of 
Mg I, Ti I, Cr I, Ni I, and to a lesser degree V I. 
The abundances derived in this way are therefore 
not direct “Fe abundances”, but an overall “metal¬ 
licity” defined by a mixture of both Fe-peak and 
a elements; we refer to this abundance as sim¬ 
ply metallicity and denote it in the remainder of 
this paper as [Fe/H], Although significantly non¬ 
solar abundance ratios, relative to Fe, could create 
additional noise in the metallicity estimates, the 
vast majority of stars in the sample are well above 
[Fe/H] ~ -2. Within this metallicity regime only 
Mg/Fe, and to a lesser extent Ti/Fe and Cr/Fe, 
are mildly non-solar and typically enhanced. This 
may introduce a small trend of slightly larger de¬ 
rived overall metallicities as the Fe abundance de¬ 
creases. 

To match the wavelength intervals defined 
above, spectra of program stars were corrected to 
zero radial velocity incorporating the RVs defined 
in § 3.1. The continuum level in the fragments of 
spectra was fit by a linear function using the con¬ 
tinuum task from the IRAF package. The pieces 
of synthetic spectra utilized for the comparison 
were convolved with a Gaussian corresponding to 
the FWHM of each echelle order to reduce their 
spectral resolution to that in the observed data. 
Then we find the best-fit model for each piece by 
the least-square method. The values and errors of 
log g and [Fe/H] were finally obtained by averag¬ 
ing these same derived parameters across the set 
of three modeled pieces. 

With this method, we derived the stellar pa¬ 
rameters (T e ff, log g , [Fe/H]) for the 434 north- 
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ern GGSS stars (see Table 1). The [Fe/H] for the 
whole GGSS sample has also been estimated pho¬ 
tometrically (Rhee et al. 2001) 3 and one can com¬ 
pare the values derived by each method (photo¬ 
metric and spectroscopic). Figure 3 shows that 
the two [Fe/H] estimates for most stars of inter¬ 
mediate metallicity ([Fe/H] = -1.0 to 0.0) follows 
a linear relation with dispersion of the order of 0.4 
dex, but systematically offset by 0.16 dex (in the 
sense that the spectroscopic metallicities tend to 
be larger). This small offset could be due partially 
to mildly elevated Mg/Fe, Ti/Fe, or Cr/Fe ratios 
in the metal poor stars. There are some extreme 
outliers in Figure 3. Some of these stars have large 
error bars, either in their spectroscopic or photo¬ 
metric metallicities. Many of the outlying points 
in Figure 3 correspond to the hottest stars in our 
sample. These stars occupy the border region be¬ 
tween giants and dwarfs on the giant-dwarf dis¬ 
crimination diagram from Majewski et al. (2000), 
and their stellar parameters determined from pho¬ 
tometric data may be biased. 

Figure 4 shows the relation between the surface 
gravity and effective temperature for the northern 
hemisphere sample. Most of the stars in Figure 
4 span ranges typical for red giants. Here we see 
the success rate of the GGSS for photometrically 
identifying bona fide red giants to be extremely 
high: more than 98 % of the stars that were photo¬ 
metrically identified by the Majewski et al. (2000) 
method to be giant stars have spectroscopic grav¬ 
ities supporting this characterization. 

3.3. Errors in Atmospheric Parameters 

In order to analyze the accuracy of the derived 
stellar parameters as a function of the quality of 
observed spectra, we performed Monte Carlo sim¬ 
ulations deriving the basic stellar parameters for 
a set of selected synthetic spectra deteriorated by 
varying amounts of added noise. We considered 
a set of signal-to-noise ratios from 3 to 30 which 
covers the typical range of S/N in our observed 
spectra. Besides random noise, we also added 
variations of T e // and RV into the model spec¬ 
tra, distributing them uniformly within ±100 K 
and ±100 m s _1 range from the true value. Fi- 

3 These photometrically derived metallicities are those that 
have been delivered by the GGSS collaboration to the SIM 
project. 


nally, we assumed that the microturbulence veloc¬ 
ity might take an arbitrary value in the range 1 to 
2 km s -1 . The stellar parameters were defined for 
a hundred resulting synthetic spectra over the typ¬ 
ical (see Figure 4) values of T e ff and log g: (4000 
K and 1.4 dex), (4400 K and 2.6 dex), and (5000 
K and 3.0 dex). The [Fe/H] took the values -2, 
-1, and 0 dex for each case. The resulting ranges 
of the estimated parameters (l-sz<?ma level) are 
shown by the solid curves in Figure 5 about the 
mean values (designated by the diamonds). The 
dashed lines represent the initial true parameters 
(log g or [Fe/H]) of the spectra. The figure shows 
that we obtain an accuracy in log g of the order of 
0.2 dex. A systematic difference in log g of ±0.1 
dex is found and is primarily the result of the dis¬ 
placement of the microturbulent velocity £ from 
the adopted 2 km s . The metallicity is defined 
with 0.1 dex accuracy and a possible displacement 
of all estimated values (again due to the unknown 
£) is -0.05 dex. To check if these results depend 
on the number of Monte Carlo simulations, we also 
evaluated some models based on sets of 400 dete¬ 
riorated spectra instead of 100. No obvious differ¬ 
ences can be seen compared with the case of 100 
simulations. 

We also extend the considered S/N ratio up 
to 100. High values of S/N do not improve the 
inferred parameters since the £ uncertainty in¬ 
troduces the most significant error. The next 
strongest factor introducing a systematic devia¬ 
tion between the "real” and estimated stellar pa¬ 
rameters is the RV error. It is seen that starting 
with S/N = 10 and better one can estimate sur¬ 
face gravity and metallicity from our spectra with 
quite good accuracy. Most of our spectra in the 
northern sample satisfy this criterion. Note that 
the accuracy of the RVs also depend on S/N. This 
can be a reason for more significant errors in log g 
and [Fe/H] obtained from our spectra. 

3.4. Stellar Parameters Derived from a 
Classical Spectroscopic Analysis 

Some of the higher S/N spectra from GGSS 
stars studied here are suitable for a detailed abun¬ 
dance analysis from which stellar parameters and 
iron abundances can be obtained via measure¬ 
ments of individual equivalent widths from a se¬ 
lected sample of Fe I and Fe II lines. This is a dif¬ 
ferent technique from the straightforward match- 
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ing of observed and synthetic spectra used in Sec¬ 
tion 3.2, which we will refer to as the “automated 
spectroscopic” method. When a strictly spectro¬ 
scopic method (using a line-by-line analysis) is 
adopted, the effective temperatures can be ob¬ 
tained by forcing a zero slope in the relation be¬ 
tween Fe I abundances with line excitation poten¬ 
tials; surface gravities are obtained from the agree¬ 
ment between the abundances from Fe I and Fe II 
lines. Another parameter that can be adjusted 
at the same time is the microturbulence velocity, 
which is tuned so that the Fe I abundances are 
independent of the equivalent widths. 

This method can provide a consistency check 
on the metallicities and log g derived with the au¬ 
tomated method presented in Section 3.2. To test 
the degree of consistency, a small sub-sample of 
10 target stars, with larger than average S/N val¬ 
ues, were selected as candidates for the detailed 
spectroscopic analysis. Our approach consisted of 
adopting the same effective temperatures from the 
Alonso et al. (1999) photometric calibration for 
the target stars and to use a sample of Fe I and 
Fe II lines that were tested to produce good re¬ 
sults for the Sun, as well as the well-studied red- 
giant Arcturus. The sampled Fe lines and their 
atomic parameters are listed in Table 2: for each 
transition the excitation potential and gf-value are 
listed. The equivalent-width measurements are 
presented in Table 3. The same spectrum analysis 
code MOOG was used in the abundance compu¬ 
tations based on the individual equivalent-width 
measurements. The adopted effective tempera¬ 
tures, and derived surface gravities, microturbu- 
lent velocities, and iron abundances are presented 
in Table 4. A comparison of surface gravities de¬ 
rived using the Fe I and Fe II lines with those 
obtained from the automated method finds a sys¬ 
tematic mean offset of [log g(automated) - log g(Fe 
I/Fe II)] = +0.4±0.4 clex. As the Fe II lines are 
quite sensitive to stellar surface gravity, this sys¬ 
tematic offset may suggest that the surface gravi¬ 
ties derived from our automated method are some¬ 
what overestimated. 

The derived Fe abundances from this detailed 
analysis are intercompared in Figure 6 with the 
metallicities obtained from Washington + DD051 
photometry (upper panel) and automated spec¬ 
troscopy (middle panel). The metallicity distri¬ 
butions shown in the middle and lower panels of 


the figure are in generally good agreement, with 
the same standard deviation of 0.18 clex, but an 
an offset of 0.18 dex in the average [Fe/H], Such 
an offset is within the quoted uncertainties in the 
metallicities obtained from the automated method 
(of roughly 0.2 to 0.3 dex). Moreover, because 
the automated method does not measure a true 
Fe abundance, but an overall general metallicity 
affected by a mixture of elements, a mild offset 
between the two methods is not unexpected. 

An important aspect of this additional verifica¬ 
tion of the results obtained with the automated 
method is a check on the derived surface gravities 
using Fe I and Fe II lines directly, via ionization 
equilibrium. As discussed previously, the GGSS 
targets were selected from Washington photome¬ 
try two-color diagrams to be giant stars and not 
dwarfs. The values of log g derived for this sub¬ 
sample of stars, as can be seen in Table 4, confirms 
their evolved status as red giants. Both the direct 
analysis of Fe I and Fe II ionization equilibrium, 
as well as the automated direct comparison with 
model spectra, indicate that the Washington + 
DD051 two-color diagram is an effective method 
for identifying red giants. 

4. RV Stability Versus Stellar Parameters 

We now focus on the 148 stars of the north¬ 
ern sample whose RVs were measured two or more 
times with large epoch differences (a time between 
successive observations of the same star not less 
than 4 months). From these data we estimate the 
RV dispersion which we designate as cr and utilize 
this as a measure of variability of the RV. Note, 
this is not the statistical standard deviation widely 
designated by cr; the present cr is calculated from 
multiple RV measures after weighting each by the 
inverse square of the accuracy of each RV measure 
( Wj ): RV m = E RVj ■ Wj / S Wj, and 

cr = ^E Wj • (RVj - RV m ) 2 / HWj. (1) 

Here Wj is the weight. If the accuracy of the corre¬ 
sponding RV was better than 50 m s -1 , we assume 
it equals 50 m s -1 because of the wavelength cal¬ 
ibration uncertainty (see §3.2). 

The derived values of cr are shown in Table 5. 
Figure 7 (upper panel) shows the distribution of 
the cr values (in km s -1 ). The largest number of 
sources occur in the first bin reveals the accuracy 
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of our estimates and confirms that its value is of 
the order of 50 m s -1 , as it follows from the accu¬ 
racy of the wavelength calibration (section §2.2). 

We next incorporate the repeated RV data ob¬ 
tained with CORALIE for the southern GGSS 
subsample (Arenas et al. 2002) in Figure 7 (lower 
panel) with 341 stars in this sample (see Table 6). 
Variation of the radial velocity a for the southern 
stars observed two and more times is estimated 
by Equation (1). A gaussian curve correspond¬ 
ing to 100 m s^ 1 sigma is shown in Figure 7 by 
dashed line. Both southern and northern subsam¬ 
ples exhibit similar behavior in their respective 
RV-variability distributions. 

We found no large, apparent systematic differ¬ 
ence in the stellar parameters between the RV- 
stable (low a) and unstable (high a) fractions of 
the sample. However, for stars with a < 100 nr 
s^ 1 the ranges of the stellar parameters appear to 
be narrower than for red giants with a > 100 nr 
s -1 . Thus, the fraction of RV-stable stars is low 
among very cool ( T e ff < 4000 K) as well as nretal- 
poor stars with [Fe/H] < -1.0, in agreement with 
previous conclusions by Jorissen et al. (1997) and 
Carney et al. (2003). 

The fraction of RV-stable and unstable stars 
can be estimated from the distributions of T e ff, 
log g , [Fe/H], and absolute magnitude My pre¬ 
sented in the left-hand panels of Figure 8; here the 
distributions of stars with a < 100 nr s -1 (dashed 
lines) are shown in comparison with the distribu¬ 
tions for all stars (dotted lines). An assessment of 
My is made with T e ff , log g, and an assumption of 
fixed mass for the red giants (we assume it equals 
to 0.9 Mq ; variation of the mass introduces in¬ 
significant scatter in My). The right-hand panels 
of Figure 8 show the fraction of stable stars in each 
bin of the histograms shown on the left. As seen in 
the figure, the fraction of stable stars has a peak 
in the range 4300-4700 K (secondary peaks at the 
wings of the distribution are caused by small num¬ 
ber statistics). Note that this range corresponds to 
the extinction-corrected ( J — K) 0 values from 0.59 
to 0.73 mag, and these values do not depend on 
the choice of color-temperature calibration. Fig¬ 
ure 2 shows that the bulk of the giant stars selected 
from the GGSS Washington+UI?051 photometry 
survey occupies this range of temperatures, show¬ 
ing that the GGSS happens to be optimized for 
selecting stars in the most RV-stable temperature 


range. 

A similar peak takes place at log g = 2.3 - 3.2 in 
the surface gravity distribution. We also find that 
the largest fraction of RV-stable stars happens for 
stars of slightly subsolar metallicity (-0.5 < [Fe/H] 
< -0.1), and then the fraction of stable stars ap¬ 
pears to drop precipitously near solar metallicity. 
Note, however, that the latter conclusion needs for 
an additional check since the most metal rich bin 
is only sampled by 6 stars. 

An analog of Figure 8 is drawn for the southern 
subsample in Figure 9. Here we use the ’’photo¬ 
metric” [Fe/H] derived by Rhee et al. (2001). The 
surface gravity is established by way of interpo¬ 
lation of the isochrones published by Bergbusch 
& Vandenberg (1992). The histograms for T e ff 
defined by the same method for both subsamples 
of northern and southern stars, are qualitatively 
similar in Figures 8 and 9. Though the other pa¬ 
rameters are estimated by different ways, the cor¬ 
responding histograms for [Fe/H], log g , and My 
are qualitatively rather similar. 

In order to better constrain what types of stars 
may be most RV-variable we again focus on the 
McDonald data and the stellar parameters derived 
from these spectra. In Figure 10 the values of log g 
are plotted versus the effective temperatures: this 
is a spectroscopic version of an HR-diagram. The 
top panel contains those stars with RV variability 
of less than 100 m s -1 , while the bottom panel are 
the definite RV variables having a > 100 nr s . 
Superimposed on these stellar points are isochrone 
curves from Girardi et al. (2000); two metallic¬ 
ity isochrones are shown (the solid curves have 
[Fe/H]=0.0 while the dashed curves have [Fe/H]=- 
0.7) with two ages for each metallicity (3.5 and 8.9 
Gyr). Plotted this way, the initial RV-results indi¬ 
cate that those stars exhibiting variability of less 
than ~100 m s _1 tend to be more concentrated in 
the log g — T e ff plane, near the He-core burning 
clump for near-solar metallicities. The stars tend¬ 
ing to show RV-variability fall all along either well- 
defined first-ascent RGB, or AGB, both at some¬ 
what low metallicities. Given uncertainties in both 
the models and the gravities and temperatures de¬ 
rived by us (which almost certainly also carry some 
systematic differences), the concentration of pos¬ 
sibly RV-stable red giants near T e //~4500K and 
log g ~2.5-3.0 coincides quite closely to the red gi¬ 
ant clump for stars having metallicities near solar, 



or slightly lower. Note that real log g may have 
systematically lower values than those shown in 
Figure 10 (see also §3.4). 

The possible association of quasi-RV stability 
with the red giant clump can be investigated fur¬ 
ther by including the southern sample, although 
we do not have spectroscopically derived gravities 
for these stars. For a relatively old population, 
however, the effective temperature of the clump 
does not vary much if the metallicity does not 
vary by much: this can been seen by inspecting 
the isochrones in Figure 10. The indications from 
the McDonald results are that the RV-stable red 
giants are only mildly metal poor with a relatively 
small dispersion in metallicity. For the 36 Mc¬ 
Donald RV-stable stars having T e ff=4500±200K 
(corresponding to the observed “clump” of points 
in Figure 10), the mean metallicity is [Fe/H]=- 
0.5±0.2. The combined northern and south¬ 
ern samples are plotted as effective temperature 
histograms in Figure 11, where the sample has 
been segregated using a. There is a clearly de¬ 
fined clump of potentially RV-stable stars with 
T e ff~4500K that probably corresponds to core-He 
burning red giants with mild metal deficiencies. 
The stars with a >100 nr s -1 , shown in the bot¬ 
tom panel, are significantly more dispersed in T e fi, 
with almost equal numbers over a range in effec¬ 
tive temperature; our suggestion is that these are 
primarily first-ascent giants or AGB stars, as well 
as intrinsically RV-stable stars (such as clump gi¬ 
ants) that have companions. The fraction of stars 
whose RV-variability is caused by components is 
not known from our observations. Demonstrat¬ 
ing the presence of companions will require much 
more repeated observations. 

Preliminary results from monitoring of photo¬ 
metric variability of a subsample of GGSS stars 
by Bizyaev et al. (2004) reveal no correlation be¬ 
tween RV- and brightness variations. However, 
more precise and simultaneous photometric and 
spectroscopic observations are needed to shed light 
on the nature of the RV variations in red giants. 

4.1. Strategies for Selecting SIM Astro- 
metric Grid Stars 

Our analysis gives insight into the optimal se¬ 
lection of giant stars for the SIM Astrometric Grid 
if time and resources are limited (as indeed they 
are) for ground-based vetting of large numbers of 


stars for the most RV-stable specimens. 

As can be seen in Figure 11, a higher fraction 
of stars with a < 100 m s _1 can be found with 
4300 A' < T e ff < 4700 A'. This is a simple, 
straightforward selection criterion resulting in an 
initial success rate of 50% in finding RV-stable red 
giants. As a point of reference, the entire northern 
sample contains 39% stable stars and the southern 
sample contains 32%. If high-resolution spectra 
are also used to derive surface gravities and metal- 
licities, the additional constraints of -0.5 < [Fe/H] 
< -0.1, and 2.3 < log g < 3.2 lead to a 81% 
and 40% success rate for the northern and south¬ 
ern samples, respectively, yet significantly narrows 
the sample. Note that the northern sample re¬ 
stricted simultaneously by all stellar parameters 
reveals the fraction of RV-stable stars that is close 
to findings by Frink et al. (2001). 

More attention to exploring the solar and 
higher metallicity ranges is warranted because if 
the relatively low RV-stable fraction we find at 
these solar-like metallicities is borne out it would 
have serious ramifications for the currently se¬ 
lected SIM Astrometric Grid sample, which is 
dominated by giant stars randomly selected from 
the Tycho catalogue. Spectroscopic metallicities 
measured for nearly 500 Tycho giants by J. Crane 
(in preparation) in the preferred magnitude range 
(V < 11) for the SIM Grid stars indicates a me¬ 
dian metallicity for this sample of [Fe/H] = -0.1, 
and perhaps larger (i.e., 0.0 dex) if the arguments 
for the need of a +0.1 dex correction to the local 
metallicity scale are applied (Reid 2002; Haywood 
2002). Thus, a significant fraction of the Tycho 
sample is at metallicities where a large fraction of 
stars may have troublesome atmospheric jitter. 

5. Conclusions 

We are carrying out high-resolution spectro¬ 
scopic observations of an all-sky sample of giants 
which was pre-selected based on the GGSS pro¬ 
gram of Washington photometry. The GGSS was 
intended to identify stars with the highest like¬ 
lihood of astrometric stability — namely, subso¬ 
lar metallicity red giants. The main purpose of 
the present study is to make an initial assess¬ 
ment of the RV-variability properties of GGSS 
stars as a function of their atmospheric charac¬ 
teristics. The RV variation was estimated for 489 
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stars with spectroscopy taken on both Northern 
and Southern telescopes. Basic stellar parameters 
were estimated for the northern subsample from 
high-resolution spectroscopy. 

A surprisingly high fraction of investigated gi¬ 
ant stars have unstable radial velocities at 100 m 
s -1 level: about 2/3 of our sample. Both of the 
samples, northern and southern, although having 
been observed independently and with different 
instruments and techniques, show similar distri¬ 
butions of RV variability. Although a number of 
obvious spectroscopic binaries are included in this 
sample, much of the low-amplitude RV variabil¬ 
ity is probably due to atmospheric motions and 
temperature inhomogeneities. 

A higher fraction of stars with a < 100 m s -1 
can be found among the objects with 4300 AT < 
T e ff < 4700 A", which corresponds to (J — K) 0 = 
0.59 0.73. If we incorporate spectroscopic metal- 
licity and surface gravity information as well, and 
select only stars with -0.5 < [Fe/H] < -0.1 and 
2.3 < log g < 3.2, it will help to identify RV- 
stable candidates more effectively, but narrows the 
sample of potential candidates substantially. From 
the point of view of minimizing expensive obser¬ 
vations, the optimal way of preliminarily select¬ 
ing RV-stable stars is to use the effective temper¬ 
ature, or more exactly, the calibration indepen¬ 
dent range of NIR colors as 0.59 < (J — AT) 0 < 
0.73. This range encompasses 44% of stars with 
a < 100ms -1 in both our southern and northern 
samples; 38.5% of the initial GGSS sample of can¬ 
didate Astrometric Grid stars (1710 of 4440) lies 
in this range. The result obtained here will aid 
in the continuing efforts to define the astrometric 
reference grid for SIM. 

The GGSS follow-up observations presented 
here have been supported by the JPL and NASA 
via grant 99-04-OSS-058. The photometric por¬ 
tion of the GGSS was funded by NAS A/JPL 
grants 1201670 and 1222563 to SRM and RJP, and 
with substantial contributions of additional fund¬ 
ing from a David and Lucile Packard Foundation 
Fellowship to SRM. VVS acknowledges support 
for part of this work from NASA via grant NAG5- 
13175. We are grateful to the staff of McDonald 
observatory for their help and support of both 
photometric and high resolution spectroscopic ob¬ 
servations, and from the Carnegie Observatories 


for similarly generous contributions for photomet¬ 
ric and medium resolution spectroscopic observa¬ 
tions of the GGSS. This publication makes use 
of data products from the Two Micron All Sky 
Survey, which is a joint project of the University 
of Massachusetts and the IPAC/Calteclr, funded 
by the NASA and NSF. 
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Table 1 

Stellar Parameters for the entire McDonald Sample of 434 GGSS stars. 


Name 

RA (J2000.0) 
(HH:MM:SS) 

DEC (J2000.0) 
(DD:MM:SS) 

T eff 

K 

[Fe/HU 

(dex) 

Dp/i 

(Kpc) 

log g sp 
(dex) 

A log g sp 
(dex) 

[Fe/H] S p 

(dex) 

A[Fe/H] sp 

(dex) 

G2358+00.31 

00:00:51.96 

00:22:16.8 

4290 

-0.3 

1.3 

4.87 

0.19 

-2.80 

0.05 

G2358+00.92 

00:01:36.38 

00:14:09.9 

4399 

-0.9 

1.1 

2.24 

0.13 

-0.42 

0.05 

G0001+00.94 

00:04:20.08 

00:27:12.2 

4324 

-0.7 

3.0 

2.67 

0.41 

-0.40 

0.16 

G0011+05.87 

00:14:18.85 

05:57:37.5 

4440 

-0.5 

0.7 

2.13 

0.19 

-0.47 

0.09 

G0011+16.75 

00:14:26.65 

17:16:03.0 

4594 

-0.8 

1.8 

1.47 

0.38 

-0.67 

0.19 


Note.—N ame of star, its RA and DEC, effective temperature, photometric metallicity [Fe/H], photometric distance, spectroscopic surface 
gravity log g and its accuracy, and spectroscopic metallicity and its accuracy. 

The full version of the Table is presented in electronic form. 
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Table 2 

The Fe I & Fe II Lines 


A(l) 

Species 

x(eV) 

log gf 

5307.361 

Fe I 

1.608 

-2.970 

5322.041 

Fe I 

2.279 

-2.840 

5497.516 

Fe I 

1.011 

-2.840 

5501.464 

Fe I 

0.958 

-2.957 

5522.447 

Fe I 

4.209 

-1.400 

5536.583 

Fe I 

2.832 

-3.812 

5539.284 

Fe I 

3.642 

-2.660 

5549.948 

Fe I 

3.695 

-2.904 

5559.638 

Fe I 

4.988 

-1.829 

5560.207 

Fe I 

4.435 

-1.040 

5577.031 

Fe I 

5.033 

-1.551 

5579.335 

Fe I 

4.231 

-2.406 

5607.664 

Fe I 

4.154 

-2.258 

5608.974 

Fe I 

4.209 

-2.402 

5611.361 

Fe I 

3.635 

-2.993 

5618.631 

Fe I 

4.209 

-1.260 

5619.224 

Fe I 

3.695 

-3.256 

5636.696 

Fe I 

3.640 

-2.608 

5638.262 

Fe I 

4.220 

-0.720 

5661.012 

Fe I 

4.580 

-2.432 

5677.684 

Fe I 

4.103 

-2.694 

5679.025 

Fe I 

4.652 

-0.770 

5696.103 

Fe I 

4.549 

-1.997 

5698.023 

Fe I 

3.640 

-2.689 

5705.466 

Fe I 

4.301 

-1.360 

5717.835 

Fe I 

4.284 

-0.980 

5724.454 

Fe I 

4.284 

-2.627 

5759.261 

Fe I 

4.652 

-2.073 

5760.345 

Fe I 

3.642 

-2.490 

5784.657 

Fe I 

3.397 

-2.673 

5793.913 

Fe I 

4.220 

-1.697 

5806.717 

Fe I 

4.608 

-0.900 

5807.782 

Fe I 

3.292 

-3.404 

5809.217 

Fe I 

3.884 

-1.690 

5811.917 

Fe I 

4.143 

-2.427 

5814.805 

Fe I 

4.283 

-1.820 

5837.700 

Fe I 

4.294 

-2.337 

5838.370 

Fe I 

3.943 

-2.337 

5844.917 

Fe I 

4.154 

-2.940 

5845.266 

Fe I 

5.033 

-1.818 

5849.682 

Fe I 

3.695 

-2.993 
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Table 2 —Continued 


A(l) 

Species 

x(eV) 

log gf 

5853.149 

Fe I 

1.485 

-5.268 

5856.083 

Fe I 

4.294 

-1.640 

5861.107 

Fe I 

4.283 

-2.452 

5000.743 

Fe II 

2.778 

-4.745 

5018.440 

Fe II 

2.891 

-1.213 

5132.669 

Fe II 

2.807 

-4.000 

5234.625 

Fe II 

3.221 

-2.240 

5264.812 

Fe II 

3.230 

-3.188 

5284.098 

Fe II 

2.891 

-3.010 

5325.559 

Fe II 

3.221 

-3.170 

5414.046 

Fe II 

3.221 

-3.620 

5425.247 

Fe II 

3.199 

-3.210 

5534.847 

Fe II 

3.244 

-2.770 


Note.- Adopted excitation poten¬ 
tials and gf-values for Fe I and Fe II 
lines 
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Table 3 

Equivalent-Width Measurements 


A #1 #2 #3 #4 #5 #6 #7 #8 #9 #10 


5307.361 

132 

109 

185 

117 

5322.041 

92 

84 

140 

86 

5497.516 


183 


172 

5501.464 




166 

5522.447 


50 


47 

5536.583 

35 

25 


17 

5539.284 

42 


54 

23 

5549.948 

32 

15 

44 


5559.638 

18 


19 


5560.207 

63 



57 

5577.031 



22 


5579.335 


14 

34 

10 

5607.664 



44 


5608.974 

25 


35 


5611.361 

31 


37 

14 

5618.631 

70 

58 

91 

56 

5619.224 

21 

11 

195 

6 

5636.696 

40 


65 

24 

5638.262 

98 

78 

126 

78 

5661.012 



15 


5677.684 

20 


24 


5679.025 

66 

54 


53 

5696.103 

27 

17 



5698.023 


28 

50 


5705.466 

50 

48 

88 

44 

5717.835 

76 

69 

105 

66 

5724.454 

19 


19 


5759.261 



20 


5760.345 

44 


65 

29 

5784.657 

52 

36 

77 

34 

5793.913 

51 

38 

71 

39 

5806.717 

60 

55 

90 

49 

5807.782 

27 


41 

16 

5809.217 

70 

54 

98 

53 

5811.917 


15 

35 


5814.805 

36 

34 

59 


5837.700 


13 


12 

5838.370 


25 


21 

5844.917 





5845.266 





5849.682 


14 

31 



118 

154 

154 

165 

116 

135 

85 

114 

100 

118 

819 

102 



73 






34 





32 

45 

45 

44 

28 





31 

21 




15 

11 



60 

72 

68 

68 

56 

62 





15 




25 

23 

16 


25 


35 

33 


27 


24 

31 



13 

19 


37 

27 


19 

62 

78 

71 

76 

58 

64 

10 

15 

21 

15 




48 


40 

32 


89 

99 

92 

98 

77 

96 



13 




15 

17 

24 

20 



63 

74 

72 

74 

51 

66 


22 



18 

18 

31 

39 

47 



30 

54 

63 

65 


47 

52 

67 

86 

81 

82 


72 





14 

9 

13 


50 

15 

12 

11 

41 

53 

23 


38 

36 

42 

54 

57 

60 

44 

41 

46 

55 

54 

57 

42 


56 

68 

66 

66 


64 

23 

26 

39 

30 

17 


56 

76 

70 

79 


70 



29 

25 

16 


37 


48 

42 


33 

16 


28 

21 



32 


45 



31 



16 


8 

8 


10 

14 


9 

10 

17 

20 

30 

26 

19 

14 


15 












Table 3 —Continued, 


A 

#1 

#2 

#3 

#4 

#5 

#6 

#7 

#8 

#9 

#10 

5853.149 


32 

67 

22 

32 

55 

58 

57 


27 

5856.083 

51 


70 


42 

51 

53 

56 

39 


5861.107 







22.7 




5000.743 











5018.440 


175 

201 



222 


186 

174 

161 

5132.669 

39 


37 

41 

35 

49 

45 



25 

5234.625 

87 

71 

100 

92 

72 

100 

86 

81 

72 

786 

5264.812 

45 

43 

55 

54 

44 

57 

56 

44 

38 

42 

5284.098 

64 

63 

82 

76 




67 

62 


5325.559 


42 

54 

56 

47 

57 

56 

45 


36 

5414.046 

30 

30 

30 

37 

34 


42 

30 

29 

20 

5425.247 

46 

41 

52 

57 

46 

61 

51 

44 

38 

38 

5534.847 


60 

70 

73 

59 

74 

70 





Note.— Equivalent widths of Fe I and Fe II lines (given in mA) for 10 stars. Stellar 
identifications are provided in Table 4. 


Table 4 

Stellar Parameters and Metallicities derived from classical spectroscopic analysis for 

10 SELECTED STARS 


Star 

ID 

T e ff 

(K) 

log g 
(dex) 

e 

(km/s) 

[Fe/H] 

(dex) 

G0858+00.192 

#1 

4630 

2.2 

1.1 

-0.31 ±0.07 

G0901+00.175 

#2 

4760 

2.6 

0.9 

-0.48 ±0.06 

G1007+00.5 

#3 

4490 

2.5 

1.8 

-0.12 ±0.05 

G1018-05.46 

#4 

4670 

1.6 

1.2 

-0.77 ±0.05 

G1030+00.15 

#5 

4780 

2.6 

0.9 

-0.34 ±0.06 

G1042-05.69 

#6 

4370 

1.1 

1.4 

-0.63 ±0.06 

G1053+00.4 

#7 

4400 

1.2 

0.9 

-0.35 ±0.07 

G1113+00.20 

#8 

4240 

1.6 

1.4 

-0.56 ±0.06 

G1147-05.31 

#9 

4600 

2.1 

0.6 

-0.48 ±0.07 

G1147-05.49 

#10 

4720 

3.0 

1.3 

-0.41 ±0.08 


Note.— Name of star, its ID in Table 3, effective tempera¬ 
ture, surface gravity, microturbulent velocity, and [Fe/H]. 
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Table 5 

RV-variability for the 148 McDonald Repeated Stars 


Name 

M 

(mag) 

T e ff 

(K) 

Dph 

(Kpc) 

[Fe/H]p h 

(dex) 

log g sp 
(dex) 

A log g sp 
(dex) 

[Fe/H] sp 

(dex) 

A [Fe/H] sp 
(dex) 

a 

(km/s) 

N 

G0118-05.38 

11.97 

4814 

0.9 

-1.0 

3.13 

0.13 

0.00 

0.05 

0.165 

2 

G0131+00.86 

11.55 

4707 

0.9 

-1.2 

2.34 

0.17 

-0.80 

0.14 

0.919 

2 

G0142+05.56 

10.91 

3959 

2.2 

-0.9 

1.20 

0.65 

-0.60 

0.28 

0.489 

2 

G0142+05.70 

12.27 

4196 

2.3 

-0.6 

1.93 

0.25 

-0.67 

0.25 

0.287 

2 

G0151+00.72 

11.29 

4334 

1.7 

-0.6 

2.34 

0.06 

-0.07 

0.22 

0.287 

3 


Note.—N ame of star, Washington M magnitude (which is approximately Johnson’s V), effective temperature, photometric 
distance, photometric metallicity [Fe/H], spectroscopic surface gravity log g and its accuracy, spectroscopic metallicity [Fe/H] 
and its accuracy, RV variability, and number of RV-observations for the northern sample objects. 

The full version of the Table is presented in electronic form. 


Table 6 


RV-variability for the 341 


CORALIE Repeat Observation Stars 


Name 

M 

(mag) 

T eff 

(K) 

D ph 

(Kpc) 

[Fe/H] p h 

(dex) 

a 

(km/s) 

N 

G0000-56.81 

11.92 

4066 

3.5 

-1.1 

0.004 

3 

G0012-28.38 

11.71 

3891 

3.8 

-1.0 

0.182 

3 

G0016-39.1207 

11.46 

4418 

4.3 

-2.0 

1.104 

7 

G0016-39.2075 

11.93 

4580 

1.7 

-1.0 

0.635 

2 

G0016-39.3290 

11.41 

4753 

0.8 

-0.7 

2.060 

4 


Note. —Name of star, Washington M magnitude, effective tempera¬ 
ture, photometric distance, photometric metallicity [Fe/H], RV variabil¬ 
ity, and number of RV observations for objects of the southern sample. 

The full version of the Table is presented in electronic form. 
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Fig. 1.— An example of the order-to-order 
radial velocity differences for the program star 
G1113+00.20. The wavelength corresponds to the 
middle of each respective order. 




[Fe/H] 

PHOTOMETRIC 


Fig. 3.— The comparison of metallicities es¬ 
timated from high-resolution spectroscopy (see 
text) and “photometric” [Fe/H] (Rhee et al. 2001). 
The dashed line shows the “robust” least absolute 
deviation fit to the data. The values of [Fe/H] 
for most of the stars with intermediate metallicity 
approximately follow a linear relation with disper¬ 
sion of order 0.4 dex. The spectroscopic [Fe/H] is 
systematically higher than the photometric one by 
about 0.16 dex. 


Fig. 2.— Distribution of T e // for the general sam¬ 
ple of 4440 GGSS candidates. 
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Fig. 4.— The relation between surface gravity and 
effective temperature for 434 stars from our north¬ 
ern sample. 




Fig. 5.— The 1-sigma level scatter of estimated 
stellar atmospheric parameters (solid curves ) as 
a function of S/N. The mean values of log g 
and [Fe/H] found for one hundred deteriorated 
model spectra are designated by the diamonds. 
The dashed lines represents the initial (true) pa¬ 
rameters (log g or [Fe/H]) of the model spec¬ 
tra. The results of simulations are shown for 
models typical of the GGSS red giants: T e ff = 
4000, log g = 1.4, T e ff = 4400, log g = 2.6, and 
T e ff = 5000, log q = 3.0, and for metallicities 
[Fe/H]=0, -1, and-2. 
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[Fe/H] classic 



- 1.2 - 1.0 - 0.8 - 0.6 - 0.4 - 0.2 
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Fig. 6. Metallicities for 10 GGSS stars measured 
three different ways: using Washington photome¬ 
try, the automated method, and a detailed spec¬ 
troscopic analysis. The three methods show gen¬ 
eral agreement. 



Fig. 7. Distribution of the variation of radial 
velocity, a (in km s -1 ) for the northern (upper 
panel) and southern (lower panel) subsamples (148 
and 341 stars, respectively) with two or more es¬ 
timates of RV. The dashed line shows a gaussian 
distribution with 100 m s^ 1 sigma). 



Mabs(v) Mabs(v) 

Fig. 8.— Left panels: histograms of distribution 
of T e ff, log g , [Fe/H], and My for the stars whose 
<j < 100 m s -1 (dashed line) in comparison with 
the whole sample (dotted line). The right side 
panels show the corresponding fraction of stable 
stars in bins of the histogram. 
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Fig. 9.— The same as in Figure 9, but for the 
southern subsample. 



Fig. 10.— The spectroscopically derived surface 
gravity, log g, plotted versus T e s for the McDon¬ 
ald northern sample: The top panel shows stars 
with a < 100 m s -1 while the bottom panel shows 
those with RV-variability greater than 100 m s -1 . 
The RV-stable stars are more concentrated in the 
log g - T e ff plane than the variables. The con¬ 
tinuous curves are isochrones from Girardi et al. 
(2000), with the dashed curves having [Fe/H]=- 
0.7 and the solid curves [Fe/H]=0.0. The two 
isochrones for each metallicity have different ages 
of 3.5 Gyr and 8.9 Gyr, respectively. The concen¬ 
trated “clump” of RV-stable stars may correspond 
to the core-He burning red giant phase of stellar 
evolution. 
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Fig. 11.— Effective temperature histograms seg¬ 
regated as the RV-stable stars (top panel) and 
RV-variable stars (bottom panel) for the combined 
northern and southern data sets. The RV-stable 
stars are sharply peaked near T e ff= 4500±200K, 
while the RV-variable stars are more spread out 
in temperature. This provides additional support 
for the suggestion that the core-He burning stars, 
with only mild metal deficiencies, may be intrinsi¬ 
cally more stable in radial velocity. 
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